Local absorption spectra of artificial atoms and molecules 
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We investigate theoretically the spatial dependence of the linear absorption spectra of single and 
coupled semiconductor quantum dots, where the strong three-dimensional quantum confinement 
leads to an overall enhancement of Coulomb interaction and, in turn, to a pronounced renormal- 
ization of the excitonic properties. We show that — because of such Coulomb correlations and the 
spatial interference of the exciton wavefunctions — unexpected spectral features appear whose in- 
tensity depends on spatial resolution in a highly non-monotonic way when the spatial resolution is 
comparable with the excitonic Bohr radius. We finally discuss how the optical near-field properties 
of double quantum dots are affected by their coupling. 
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I. INTRODUCTION 

In recent years much attention is being devoted to 
the properties of semiconductor quantum dots (QDs). 
In these systems, carriers are subject to a confining po- 
tential in all spatial directions, giving rise to a discrete 
energy spectrum ("artificial atoms") and to novel phe- 
nomena of interest for fundamental physics as well as 
for applications to electronic and optoelectronic devices 
iQ,^ . The extension and the shape of the QD confining 
potential varies, depending on the nanostructure fabri- 
cation technique: The dots that are studied most ex- 
tensively by optical methods are induced by quantum 
well (QW) thickness fluctuations or obtained by 

spontaneous island formation in strained layer epit axy 
10-^ , self-organized growth on patterned substrates , 
stressor- induced QW potential modulation jl^, cleaved 



edge overgrowth 
techniques p3|,p| 



|13], as well as chemical self-aggregation 
The resulting confinement lengths fall 
in a wide range between 1/im and 10 nm. 

In spite of the continuing progress, all the available 
fabrication approaches still suffer from the effects of in- 
homogeneity and dispersion in the dot size, which lead to 
large linewidths when optical experiments are performed 
on large QD ensembles. A major advancement in the 
field has come from different types of local optical experi- 
ments, that allow the investigation of individual quantum 
dots thus avoiding inhomogeneous broadening [p|-p^ . 

Among local spectroscopies, the approaches based on 
scanning near-field optical microscopy (SNOM) ||l^ are 
especially interesting as they bring the spatial resolu- 
tion well below the diffraction limit of light: With the 
development of small- aperture optical fiber probes, sub- 
wavelength resolutions were achieved (A/8 — A/5 [|l^ or 
A/40 Q) and the first applications to nanostructures 
became possible 18-2^. As the resolution increases, 
local optical techniques in principle allow direct access 
to the space and energy distribution of quantum states 
within the dot. This opens, however, a number of ques- 
tions regarding the interpretation of these experiments, 



that were often neglected in the past. 

First of all, for spatially inhomogeneous electromag- 
netic (EM) fields it is no longer possible to define and 
measure an absorption coefficient that locally relates the 
absorbed power density with the light intensity (since 
the susceptibility x(r, r') cannot be approximated by a 
local tensor). In the linear regime, a local absorption 
coefficient can still be defined, which is however a com- 
plicated function that depends on the specific EM field 
distribution The interpretation of near- field spec- 

tra therefore requires calculations based on a reasonable 
assumption for the profile of the EM field. 

Secondly, the quantum states that are actually probed 
are few-particle states of the interacting electrons and 
holes photoexcited in the dot. Even in the linear regime, 
excitonic effects are known to dominate the optical spec- 
tra of dots since Coulomb interactions are strongly en- 
hanced by the three-dimensional confinement. Near-field 
spectra probe exciton wavefunctions, and their spatial 
coherence and overlap with the em-field profile will de- 
termine the local absorption p^ ]. 

In this paper, we show how the above phenomena affect 
local spectra of QDs, with special attention to the case 
of coupled dots ( "artificial molecules" ) where carriers in- 
teract across the barrier via tunneling and/or Coulomb 
matrix elements Indeed, the optical properties of 

coupled dots are currently of great interest not only in 
view of the unavoidable inter-dot interactions occurring 
in real samples with dense QD packing, but also in view 
of their relevance for designing novel devices including 
gates for possible solid-state implementations of quan- 
tum information processing [p6| . 

We will show that the relative phase of the exciton 
wavefunction in adjacent coupled dots (or in different re- 
gions of the same dot) can induce dramatic changes in 
the selection rules with respect to far-field spectra: A 
realistic prediction of these effects require accurate cal- 
culations taking into account quantum confinement as 
well as Coulomb interactions. Our theoretical scheme 
is especially designed to allow a realistic description of 
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the quantum states of the interacting electron and holes 
photoexcited in the linear regime. In this respect we 
improve drastically over previous approaches, which gen- 
erally focused on a more detailed treatment of the EM 
field distributions 

Our theoretical framew ork for dots is summarized in 
Section ||, while Sections |in| and 



IV discuss our results 



and conclusions for single and coupled dots. 



II. THEORY 



In this section, we summarize our theoretical approach 
for computing local absorption spectra for semiconduc- 
QDs. 



tor QDs. We first show in Sec. [I A how to compute 
the single-particle eigenstates for electrons and holes 
subjected to a three-dimensional confinement p otent ial. 
These single-particle states are then used in Sec. [IB for 
the calculation of electron-hole (i.e., optical) excitations. 
In analogy to semiconductor systems of higher dimen- 
sionality, we shall refer to these excitations as excitons; 
the properties of such excitons, however, are not only 
governed by the attractive electron-hole Coulomb inter- 
action but in addition by the strong quantum confine- 
ment. Finally, we use in Sec. [I C the above ingredients 
to derive the equations needed for the calculation of local 
optical absorption spectra. 



B. Exciton states 

When the dot structure is perturbed by an external 
light field (e.g., laser), electron- hole pairs are created 
which propagate in the presence of the mutual Coulomb 
interaction and of the dot confinement potential. Within 
the present paper, we shall restrict ourselves to the linear 
optical response, i.e., the dynamics of a single electron- 
hole pair. Then, the exciton dynamics is described by the 
electron- hole wavefunction '^{r^jTh), with the squared 
modulus being the probability of finding the electron at 
position Tg when the hole is at position Vh- 

If we expand the electron- hole ( "exciton" ) eigenfunc- 
tion in terms of single-particle states viz.: 

vl/^(re,r,) = ^</.^(r,)vl/;:>::(r,0, (2.2) 
we obtain the excitonic eigenvalue problem [B5,p5[ : 



g'')^^ +\^V'''\ (2 3) 

^' v' 

As will be shown in the following, the exciton spec- 
trum E\ directly provides the optical transition energies 
whereas the excitonic wavefunctions determine the 
oscillator strengths of the corresponding transitions. In 
Eq. (2.3) we have introduced the electron-hole Coulomb 
matrix elements Ip6| : 



A. Single-particle states 



In semiconductor QDs, carriers are confined in all 
three space directions. To simplify our analysis, we as- 
sume that a suitable parameterization of the dot con- 
finement potential is known (e.g., from experiment) and 
that the confinement potential varies sufficiently slow on 
the length scale of the lattice constant. We thus shall 
make use of the envelope- function approach [3^ ]; more- 
over, since the energy region of our present concern is 
relatively close to the semiconductor band gap, we de- 
scribe the material band structure in terms of a single 
electron and hole band within the usual effective-mass 
approximation. More specifically, the envelope-function 
equation for single electrons and holes reads: 



r_|!Z! + v:'\v)] r,'\r) = 6f (r), (2.1) 

where irie {mh) is the effective mass and {^1^) is the 
confinement potential energy for electrons (holes). Fol- 
lowing our app roac h developed earlier we numeri- 
cally solve Eq. (|]^) for arbitrary confinement potentials 
by use of a plane- wave expansion with periodic boundary 
conditions (see Appendix A). 



v: 



eh 



dr^dr 



h 



0^*(re)</>^,(re)0r(r;,)e(r/.) 



(2.4) 



where e is the elementary charge and is the static di- 
electric con stant of the bulk semiconductor (note that 
in Eq. (2.4) we have not considered the electron-hole 
exchange interaction). Within our computational ap- 



proach, we consider in Eq. (2.3) typically a basis of 12 



states for electrons and holes, respectively, and obtain 
the excitonic eigenf unctions by direct diagonalization of 
the Hamiltonian matrix. 



C. Local optical absorption 

When the semiconductor nanostructure is excited by a 
local near-field probe, the total absorbed power a{uj) at 
a given frequency lu is proportional to J dr £^{r)P{r, uj), 
where £uj{tc) is the electro- magnetic field distribution of 
the near-field probe. Within linear response, the in- 
duced interband polarization P{r,uj) is related to £'aj(r) 
through: 



P(r,c^) = J dr'x(r,r';c^)£^(r'), 



(2.5) 
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where the non-local electrical susceptibility x(r, r'; can 
be expressed in terms of the excitonic eigenenergies and 
eigenf unctions M: 



X(r,r';w) 



A 



\]/^(r,r)\E'^*(r',r') 



E\ — flLO 
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(2.6) 



Here, is the dipole-matrix element of the bulk semi- 
conductor, while we have introduced a small damping 
constant 7 accounting for the finite lifetime of exciton 
states due to environment coupling (e.g., phonons). To 
derive our final expression, it turns out to be convenient 
to consider for the elctromagnetic field distribution a 
given profile ^ centered around the beam position R, i.e., 
Sui{y) — £^^(r — R). Then, the local spectrum for a given 
tip position R can be expressed in the form ||2^ (sec also 
Appendix A): 



Ex — — huj 



(2.7) 



where 



a^(R) 



dr «'^(r,r)^(r-R) 



(2.8) 



Two limiting cases can be identified. For a spatially ho- 
mogeneous electromagnetic field (far- field), the oscilla- 
tor strength is given by the spatial average of the 
excitonic wavefunction, i.e., a^(R) — \ J dr \E''*'(r, r)p. 
In the opposite (and hypothetical) limit of an infinitely 
narrow probe, ^(r — R) = 6{r — R), one is probing the 
local value of the exciton wavefunction, i.e., a^(R) = 

|\E''^(R, R)p. Finally, within the intermediate regime of 
a narrow but finite probe, ^'^(r, r) is averaged over a re- 
gion which is determined by the spatial extension of the 
light beam; therefore, excitonic transitions which are op- 
tically forbidden in the far-field may become visible in 
the near-field. 



III. RESULTS 

In the following sections we consider the interaction of 
the EM field with excitonic states of single and double 
QDs; for the latter system, we focus particularly on the 
transition between two isolated QDs and a the "artificial 
molecule", where the electronic states of two QDs are 
strongly overlapping. 



A. Single-particle states 

We shall consider a prototypical QD confinement 
which is composed of a 2D harmonic potential in the 



(x, ?/)-plane and a rectangular quantum well along z; such 
confinement potentials have been demonstrated to be 
a good approximation for self-assembled quantum dots 
formed by strained-layer epitaxy. We focus on cases 
where the z-confinement is stronger than the (x, y) one, 
so that the confinement potential can be written as 



^'•\x,y,z) = V^^'\x,y) + V,''H{\z\ - |), (3.1) 

where Zq is the width of the quantum-well and VJ^^'' the 
band offsets for electrons and holes, respectively. For a 
single dot the in-plane confinement potential V^''^{x^y) 
is of the form: 



V;-'\x,y) = -K,,,,{x^+y^), 



(3.2) 



while for two dots (i.e., double dot) separated by the dis- 
tance d: 



y\^'''{x,y) = 



lK.e.h{{\x\-if + y^) ior\x\> 
i/Ce,,((f -a:2)+y2) 



X 

otherwise 



(3.3) 



with JCe.h — mejii'^Jt'^Y' 1 and fiLo'p'^ the level splittings 
of the in-plane harmonic potential. The shape of the 
double-dot potential has been obtained by matching the 
parabolas with opposite curvature, such that the poten- 
tial is continuous and smooth at a; = the shape of the 
resulting potential along the x- direction is shown, for se- 
lected inter-dot distances d, in Figs. 1(b), 1(d). Material 
and dot parameters which are used in this paper are listed 
in Table |; with this choice of parameters, electron and 
hole wavcfunctions have approximately the same lateral 
extension, and the QW-induced intersubband splittings 
are much larger than huo'^ and fiLOo- 

With our choice of the confinement potential, 
Eq. (refeq:total. confinement), the single-particle energies 
of a QD are Eqd = -Eqw + -Eha,m , where Eq^i is the con- 
finement energy of the QW along z and i?ha,m is the con- 
fienment energy of the 2D single- or double-harmonic po- 
tential. Single-particle energies and envelope functions 
have been computed numerically within a plane-wave 
scheme. However, for a single QD the 2D eigenstates 
can be found analytically and are the well-known "Fock- 
Darwin" states |l| (we stress, however, that the exten- 
sion in the z-direction is of crucial importance for the 
calculation of the Coulomb matrix elements and the op- 
tical properties, and unavoidably has to be taken into 
account in any realistic calculation; see also discussion in 
Ref. 0). For such states E'harm = (n -h l)hLol''^, where 
n = 0, 1, . . . is the principal quantum number, and each 
level is {n + l)-fold degenerate; in Table || we summarize 
for convenience some properties of these "Fock-Darwin" 
states p. 
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Figure 1 shows the calculated single-particle energies 
for electrons and holes for the more complex case of 
a double QD with the confinement potential given in 
Eq. ( |3.3D and with parameters listed in Table ||. The lower 
panels (Figs. l(b,d)) show the confinement potentials for 
electrons and holes at selected inter-dot distances. Obvi- 
ously, for large dot separations c? (d ^ 60 nm) the system 
can be well approximated by two separate QDs; in this 
regime the equidistance of the excited states and the cor- 
rect degeneracy of the Fock-Darwin states is obtained. 
When d is small enough that carriers have sufficient en- 
ergy to overcome (or tunnel through) the barrier between 
the two dots, the degeneracy is removed, and the energy 
levels have a non-monotonous behaviour which reflects 
the transition from two separated carrier systems to a 
single one, and is similar to the one found, e.g. for cou- 
pled QWs . For the smallest dot distances the double- 
dot potential merges into a single-dot potential, and the 
Fock-Darwin states of a single dot are recovered. 



B. Role of the Coulomb correlation in the far-field 
spectra 

Before turning to the analysis of near-field spectra, we 
shortly discuss the limiting case of very broad EM field 
distribution (far-field spectra). This discussion allows us 
to elucidate the role of the electron-hole Coulomb cor- 
relation, particularly in the transition from two separate 
"artificial atoms" to an "artificial molecule" . 

F ar-fiel d spectra can be obtained in the formalism of 
Sec. II C with a spatially homogeneous electromagnetic 



external inhomogeneous EM-field (as will be discussed 
later), the selection rules noted above are relaxed. In- 
deed, when d is reduced and the two QDs begin to inter- 
act, the calculated spectra show a much richer structure, 
as shown in Fig. 2(a) and 3(a), reflecting the reduction 
of built-in symmetry. Obviously, when c? ~ 0, the usual 
selection rules of a single-QD are recovered. 

When Coulomb interaction is included, inspection of 
the exciton wavefunctions ^'^j, (obtained from the so- 



field distribution probe ^): 



lutions of Eq. ( |2.3| )) shows that a number of different 
single-particle transitions contributes to each excitonic 
state Coulomb interaction affects the optical spec- 
tra [see Figs. 3(a) and 3(b)] in several ways. Firstly, be- 
cause of the attractive electron-hole interaction leading 
to the "bound" excitonic states, we observe a red-shift 
of the peaks; peak separation is barely affected, however, 
at least for large d. Secondly, we observe a redistribu- 
tion of the oscillator strength; for isolated QDs the three 
main peaks are of similar height. In general, oscillator 
strength is transferred from higher to lower peaks. This 
effect is particularly strong, e.g., in the doublet which 
splits from the lowest peak when the two QDs approach; 
contrary to the uncorrelated case, the heighest partner 
is very weak. Finally, Coulomb interaction is responsible 
for the appearance of additional lines (see, e.g., for d — 70 
nm the peak at 70 meV) . While the first two effects (red- 
shift and transfer of oscillator strength) are similar to 
what is found in the absorption spectra of semiconduc- 
tor quantum wires, and thus can be considered as a gen- 
eral fingerprint of Coulomb correlations in the optical 
properties of semiconductor nanostructures, the origin of 
the additional peaks is best discussed in connection with 
the calculated near-field optical spectra and, therefore, is 
postponed to the next section. 



^(r)=const 



dr ^'^(r,r) 



(3.4) 



Figs. 2 and 3 show the calculated far-field spectra for 
a double QD as a function of the dot distance d. We first 
concentrate on the calculations where Coulomb correla- 
tions were artificially set to zero (Figs. 2(a) and 3(a)): 
Because of symmetry, only a small fraction of all p ossi- 
ble electron-hole transitions is visible; from Eq. (^^) and 
using / c?(/?expi(me -l- irih)^ oc dmc,-mh obtain that 
optical transitions are only allowed between electron and 
hole single-particle states with opposite angular momen- 
tum. Indeed, for large distances (uncoupled QDs) only 
three strong absorption peaks are observed, with an en- 
ergy splitting of approximately huj^ + Ticu^; the intensity 
of the peaks increases with energy (with ratio 1:2:3). 
These can be attributed to transitions between electrons 
and holes single-particle states (see Table ||) of the Is 
symmetry (peak at « 70 meV), the Ip symmetry (peak 
at « 95 meV), and the Id and 2s symmetries (peak at 
w 120 meV). 

When symmetry is reduced, either because of an asym- 
metric confinement potential or by the presence of an 



C. Optical near- field spectra 

In this section we discuss the local absorption spectra 
of single and coupled QDs. Because of the narrow well 
width of the dot confinement potential (see Table |), the 
EM profile of the near-field probe along z has an only 
minor influence on the results, and we use: 



£.{x,y,z) cx exp ^- 



2a' 



(3.5) 



The spatial resolution of the electromagnetic fleld distri- 
bution of Eq. (B.f;) is then approximately given by the 
full-width at half maximum (FWHM) of the Gaussian 
(i.e., 2\/2 ln2(T w 2.35cr). Since the Gaussian acts as 
an envelope on 5''^, in the intermediate regime of a nar- 
row but flnite a the spatial average only extends over the 
region where the Gaussian is non- vanishing. 

Since the extension of the quantum states under in- 
vestigation is of the order of a few tens of nano-meters 
(see also Figs. 5 and 8, to be discussed below), in our 
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calculations we consider three different regimes of spatial 
resolution: i) a regime where the FWHM is much larger 
than the extension of the quantum states (as a charac- 
teristic value we use cr = 50 nm); iij a regime where the 
FWHM is comparable to the extension of the relevant 
quantum states (we use cr = 10 nm); Hi) a regime with 
an extremely narrow probe beam (we use cr = 0.1 nm). 
Calculations performed in this latter (unphysical) regime 
are used for illustrative purposes to obtain a "cartogra- 
phy" of the exciton wavefunction, as discussed at the end 



of Sec. [I C. We finally notice that the excitonic Bohr ra- 
dius w 12 nm for GaAs. 



1. Single quantum dot 

In Fig. 4 we report the calculated local absorption 
spectra a^{X,hu}) for a single QD as a function of the 
tip position. The tip is swept along one direction, pass- 
ing through the center of the QD. 

In Figs. 4(a-c) we show the calculated spectra neglect- 
ing Coulomb interaction. For the highest spatial reso- 
lution [Fig. 4(a)], the local absorption at photon energy 
E\ is proportional to J dz \I''''(r, r)|j,=o- Given the en- 
ergy splitting TiLd^ = 3.5 meV for holes and Tiu'^^ = 20 
meV for electrons, we can attribute the triplet of peaks at 
~ 70 meV to the single-particle transitions involving the 
Is state of electrons and the Is, Ip, and (2s, Id) states 
in order of increasing energy) of holes (see also Table 

); analogously, the triplet at « 90 meV is attributed 
to the transitions involving the Ip state of electrons and 
the Is, Ip, and (2s, Id) hole states; indeed, in Fig. 4(a) 
the localization of the absorption peaks is suggestive of 
the S-, p- or d-type symmetry of the corresponding Fock- 
Darwin states. These features are still present at the 
intermediate resolution [Fig. 4(b)], but disappear at the 
opposite limit of a broad probe [Figs. 4(c)]. This is ex- 
pected, since, when a localized EM-field is present, the 
symmetry of the whole system (nanostructure-(-EM-field) 
is lower than that of the nanostructure (except when the 
probe is centered in the symmetry center of the struc- 
ture) , and far-field selection rules are relaxed. When the 
probe is broadened, however, the built-in symmetry of 
the structure is recovered, and optical far-field selection 
rules (i.e., optical transitions only between electron and 
hole states with opposite angular momentum to) apply; 
therefore, the spectra are almost identical to those of two 
separated dots in far-field spectroscopy, already discussed 
in Figs. 2 and 3. 

When we compare Figs. 4(a"c) with Figs. 4(d-f), we 
find that Coulomb interaction induces several effects 
which are expected on the basis of the discussion of the 
far-field spectra. In particular, we find i) a.n almost 
rigid redshift of the spectra; ii) a. transfer of oscillator 
strength from transitions at higher energies to those at 
lower energies; Hi) the appearance of new features in the 
optical spectra. To discuss the origin of these new opti- 



cal features caused by Coulomb interactions, let us con- 
sider, e.g., the optical peaks at photon energy « 65 meV 
(Figs. 4(d-f)): They are quite strong at a = 0.1 nm 
(Fig. 4(d)), almost disappear at cr = 10 nm (Fig. 4(e)), 
and are visible again in the far field limit (Fig. 4(f)). 
Such a behaviour is rather unexpected and noticeably 
differs from that of other transitions, which — with in- 
creasing cr — either remain strong or gradually disappear 
due to symmetry reasons, as discussed above. To inves- 
tigate the origin of this non-monotonic dependence, in 
the following we analyze the three excitons within the 
corresponding energy range. Figure 5 shows a contour 
plot of the respective exciton wavefunction '^^{r,r)\z=o- 
Apparently, in Fig. 5(a) the exciton has s-type symme- 
try, whereas the other two electron-hole states have p- 
type symmetry. (Because of the periodicity box used 
in our calculations, the two- fold degenerate p-type ex- 
citon wavefunctions have cartesian rather than cylinder 
symmetry; note that, since the presence of the near-field 
tip destroys the cylinder symmetry, the wavefunctions 
shown in Fig. 5 indeed form a natural basis; see also Ta- 
ble |l^. Next, we note that the average / dr ^''''(r,r) 
of the p-type exciton wavefunctions is zero. Since with 
increasing a the radius within which the exciton eigen- 
functions "if^ are averaged increases, we expect for these 
p-type functions with increasing a a monotonically de- 
creasing behavior. The exciton shown in Fig. 5(a), on 
the other hand, has a non-zero average and is therefore 
visible in both the optical far- and near-field. A closer 
inspection of the exciton wavefunction 5*^^ reveals that 
the largest contribution stems from the transition be- 
tween the Is state of electrons and the 2s state of holes, 
but there is also a noticeable contribution from the Is- 
Is and Ip-lp electron-hole transitions. Indeed, only the 
latter contributions couple in the far-field to the light 
field. In the regime of finite resolution, there is an opti- 
mal cancellation when the FWHM of the EM near field 
(i.e. 2\/2 ln2cr « 2.35cr) becomes equal to the Bohr ra- 
dius. This is clearly depicted in Fig. 6, where, in order 
to facilitate our discussion, we have introduced the quan- 
tity (X J dR (R) , which provides a measure of the 
relative contribution of each exciton to the absorption 
spectra. Fig. 6 shows for the three excitons (shown 
Fig. 5) within the energy region of 65 meV: We observe 
that with increasing cr, the p-type functions (open cir- 
cles) indeed vanish monotonically, whereas for the s-type 
exciton (full circles) there exists an optimal cancellation 
when the FWHM of the EM field distribution becomes 
approximately equal to the Bohr radius. In spite of the 
specific carrier states of a single parabolic QD, we expect 
that such non-monotonic behavior appears quite gener- 
ally in semiconductor nanostructures where carrier states 
are confined on a length scale comparable to the Bohr ra- 
dius, and thus provides a striking fingerprint of Coulomb 
correlations in the optical near-field spectra (we find sim- 
ilar behavior in our calculations for the near-field spectra 
of coupled QDs discussed below). 
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2. Double quantum dot 
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In Figure 7 we show the calculated local absorption 
spectra a^{X, Uuj) for a double QD for selected values of 
the interdot distance and a. The tip of the probe is swept 
along the direction which passes through the centers of 
the two QDs. 

Let us first concentrate on the results with a — 0.1 
nm and with the Coulomb interaction taken into account 
(Figs. 7(d,g,j)). With decreasing interdot distance we ob- 
serve the transition from a system where the energetically 
lowest exciton states are almost localized in the spatially 
separated minima of the two dots, to a system where the 
electron-hole states extend over the whole nanostructure. 
Here, the s-like ground-state excitons of Fig. 7(j) split up 
into a "bonding" and an "anti-bonding" state (Fig. 7(d)). 
By comparing Figs. 7(d) and 7(f), we find that in the op- 
tical far-field only the symmetric ground state exciton 
couples to the light field. 

Next, we discuss the optical features at the photon en- 
ergy of « 70 meV for d = 40 nm. As in the case of 
the single dot, these features show a non-monotonic de- 
pendence on the probe width. As can be inferred from 
the calculations with a = O.l nm, there are several exci- 
tonic states contributing to the spectral features in this 
energy range; Fig. 8 shows the excitonic wavefunction of 
two states out of the six states with E\ « 70 meV for for 
d = 40 nm; it can be inferred that for a spatial resolution 
of the near-field probe comparable to the excitonic Bohr 
radius (~ 12 nm) there is again an optimal cancellation. 
This is a remarkable finding, because it clearly demon- 
strates that such a behaviour indeed is a general charac- 
teristics of semiconductor nanostructures, and does not 
depend on peculiar symmetries of the confining potential. 



IV. SUMMARY AND CONCLUSIONS 



We thank Fausto Rossi for most stimulating discus- 
sions. This work was supported in part by INFM through 
PRA-99-SSQI, and by the EC under the TMR Network 
"Ultrafast Quantum Optoelectronics" and the 1ST pro- 
gramme "SQID" . U.H. acknowledges support by the EC 
through a TMR Marie Curie Grant. 



APPENDIX A: PLANE- WAVE APPROACH 

In this Appendix we discuss details of our numerical 
solution schemes based on a plane-wave expansion. Fol- 
lowing our approach developed earlier [Q, we consider 
the problem of a single or double QD which is located in- 
side a box with periodic boundary conditions, where the 
boxsize is chosen sufficiently large to avoid interactions 
with "neighbour" dots. As a complete set of functions, 
inside the periodicity box we use a plane- wave basis, |k), 
with: 



ka = ° , UaeZ, a = x,y,z. (Al) 

La 

Here La denotes the sizes of the periodicity box (we use 
the same box for electrons and holes). We next expand 
the single-particle wavefunctions for electrons and holes 
within the plane-wave basis: 

j dre-^''-0f (r), (A2) 

with 17 the volume of the periodicity box. The envelope- 
function equation (pj) is then transformed to: 



We have analyzed theoretically the interaction between 
a model near-field probe and a zero-dimensional het- 
erostructure: Quantum confinement of the electron and 
hole states, as well as their Coulomb interaction in the 
linear regime are fully included in our description. 

We have specifically considered single and coupled 
semiconductor quantum dots, and shown that absorp- 
tion is strongly influenced by the spatial interference of 
the exciton wavefunctions, which depends on the spatial 
extension of the light beam. As a consequence, near-field 
experiments on quantum dots are predicted to display 
unexpected spectral features whose dependence on spa- 
tial resolution is highly non-trivial. 

When combined with an appropriate choice of the EM 
field distribution, our approach provides the necessary 
tool for interpretation of near-field absorption spectra of 
quantum dots as the spatial resolution of experiments be- 
comes comparable with the Bohr radius of the exciton in 
the nanostructure. 



k' ' 

which can be solved by standard diagonalization tech- 
niques. To keep the numerics tractable, only wavevectors 
smaller than a given cut-off wavevector are considered 
(typically 2000-3000 wavevectors). In our computational 
approach, we perform the Fourier transform of the con- 
finement potential by storing V^''^{y) on an appropriate 
grid (with a typical number of 30 points along each di- 
rection), and approximating within each cube V^'^{y) by 
its average value. 

In the calculation of the near-field spectra, we define 
the electron- hole index I — {pL^v). Then: 

*\r,r)=^*f 0^,(rX(r), (A4) 
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and we obtain for a^(R) of Eq. (2^) the final result: 



I k,k' 

with ^k(R) = J dr ^(r)e*'^ ("-+^). 



^M.k^'.k'P, (A5) 
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TABLE I. Material parameters for GaAs/AlGaAs and dot 
parameters which were used in the calculations. Here, nio is 
the free-electron mass. 



description 


value 


units 


electron mass me 


(J.UO( 


mo 


hole mass 


0.38 


ruo 


dielectric constant Ko 


12.9 




conduction-band offset for electrons Vo 


300 


meV 


valcncc-band offset for holes 


200 


meV 


confinement energy Tuo^ for electrons 


20 


meV 


confinement energy hu)^ for holes 


3.5 


meV 


quantum-well width Zo 


10 


nm 



TABLE IL Eigenfunctions (Fock-Daxwin states) with lowest energies for a particle with mass jj, and for a potential of 
the form V{x,y) = i/LtWo(a;^ +2/^) = |/^o^^ (i-e-, two-dimensional harmonic oscillator). We use X = x/ao, y = y/ao, and 
TZ = r/tto, with Oo — \/J^^)- Because of cylindrical symmetry, the angular momentum in the ^-direction is a good quantum 
number (m) and the angular part of the wavefunctions is of the form oc exp ±invp; we use the notation s for m = 0, p for 
m = ±1, and d for m = ±2. 





cartesian coordinates: 


cylinder coordinates: 




Energy (TuVo) 


^{x,y)<xeKp-^{x^ + y^) 


0(72., (fi) oc exp — 172,^ 


notation 


1 


xl 


xl 


Is 


2 


xX 








xy 


xTZ exp ±iip 


Ip 


3 


xxy 

x{2X^ - 1) 

x{2y^ - 1) 


x{n^ - 1) 


2s 




xTZ^ exp±2i(p 


Id 
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FIG. 1. Single-particle confinement energies as a function of 
form of the confining potential along x-axis (lower panels) for d 
(dotted line). Left and right panels correspond to electrons and 



the distance between the two dots, d (upper panels) and the 
= 20 nm (solid line) , d = 30 nm (dashed line) and d = 40 nm 
holes, respectively. 
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without Goulomb interactions 
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FIG. 2. Optical absorption spectra for a honiogerieous electromagnetic field profile (i.e., far field) for a double quantum dot 
and for different distances d: (a) Coulomb interactions neglected; (b) Coulomb interactions included. We use 7 = 1 meV. The 
photon energy is measured with respect to the bandgap. 
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FIG. 3. Same as Fig. 2; the size of each dot corresponds to the height (i.e., oscillator strength) of the corresponding absorption 
peak. 
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FIG. 4. Local absorption spectra a^{X,hu}) for a single QD with (Figs. 4(d-f)) and without (Figs. 4(a-c)) Coulomb inter- 
actions and for different values of a. Photon energy hu is measured with respect to the bandgap, and X is the position of the 
tip along the x-axis (F = 0). In this calculations we use a basis of 6 electron and hole states, respectively. 
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FIG. 5. Contour plot of the cxciton wavefurictiori ^''^(r, r) for three excitons which contribute to the absorption peak at Ri 65 
meV. Sohd and dashed hnes correspond to positive and negative values, respectively. 
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FIG. 6. The relative contribution, 7j , as a function of a, for the excitons (depicted in Fig.5) which are responsible for the 
non-monotonic behavior of the feature at 65 meV. 
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FIG. 7. Local absorption spectra o:^{X, Tiuj) for a double QD with (Figs. 7(d-l)) and without (Figs. 7(a-c)) Coulomb interac- 
tions and for different values of a and interdot distance d. Photon energy huj is measured with respect to the bandgap, and X 
is the position of the tip along the x-axis {Y = Q). In our calculations we use a basis of 12 electron and hole states, respectively. 
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FIG. 8. Contour plot of the cxcitori wavcfurictiori >]/'^(r,r) of the two excitons which are responsible for the non-rnonotonic 
behaviour of the features at « 70 meV at the interdot distance d = 40 nm. Solid and dashed lines correspond to positive and 
negative values, respectively. The upper (lower) panel refers to exciton with energy 69.1 meV (70.3 meV). 



